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Abstract 

A relativistic theory of the Zeeman splitting of hyperfine levels in two-fermion systems 
is presented. The approach is based on the variational equation for bound states derived 
from quantum electrodynamics [1]. Relativistic corrections to the g- factor are obtained up 
to O ((a) 2 ). Calculations are provided for all quantum states and for arbitrary fermionic 
mass ratio. In the one-body limit our calculations reproduce the formula for the ^-factor 
(to O ((Z«) 2 )) obtained from the Dirac equation. The results will be useful for comparison 
with high-precision measurements. 

1. Introduction 

In a recent paper [2] we have presented a self-consistent variational method for calculating 
the non-relativistic Lande g factor of the two-fermion bound-state system. In the lowest- 
order approximation the linearly dependent part of the energy splitting for a two-fermion 
system placed in a weak static magnetic field B can be written as [2-6] 

AEe j%j,s,e, Sl ,s 2 = i»Bigi + HB292) Bmj, (1) 

where J, mj, S, £, Si, s 2 are quantum numbers, which characterize the system: si and s 2 are 
the spins of the first and second particle respectively, S = Si + s 2 , Sj + s 2 — 1, |«i — s 2 \ 
is the total spin of the particles, £ and J represent the orbital and total angular momentum 
quantum numbers, where J = £ + S, £ + S — 1, \£ — S\. The projection of the total angular 
momentum on the B field direction is labeled by mj = — J, — J+ 1, J — 1, J- The "Bohr 
magnetons" for the two particles are defined as ubi = Q\h/2ra\C, and fig 2 = —Q 2 h/2m 2 c, 
where Qi and Q 2 represent the magnitude of the charges. In our notation, mi and m 2 
correspond to the masses of the light and heavy particle respectively. The description of the 
interacting system by the set of quantum numbers J, mj, S, £, si, s 2 corresponds to the 
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LS coupling representation. This representation is used in contrast to the customary 
coupling scheme (for the case m 2 >> mi), where the states are taken to be the eigenstates 

of the operators ]\ = (l + s^ , j| = s|. As discussed previously [7], for the general 

case of arbitrary mass ratio, the j\ value is not a good quantum number. Even in the LS 
representation, for spin-mixed states the orbital angular momentum L, and total spin S of 
the system are not conserved. In this case we designate the states by an additional quantum 
number s, which takes on the values of or 1 for quasisinglet (sg q ) and quasitriplet (tr q ) 
states respectively. 

In our calculations we assume that the energy-level splitting (1) is smaller than the 
hyperfine structure (HFS) splitting, AE ext « AE HFS , i.e., we treat the interaction with 
an external magnetic field as a perturbation. 

The non-relativistic Lande factors g\ and g 2 obtained in [2] can be summarized as follows: 



for £= J- 1: 



for 1 = J + l: 



I'i^-A ■ i :J T ' ) "7" (2) 



j v 2 / J 

for spin-mixed states £ = J ^ 

where i — 1,2 is the index designating the particle. The index j is defined as j = 1 when 
i = 2, and j = 2 when i — 1. The quantities fa represent the mass factors: 

m i 

to = : , 5 

m 1 + m 2 

and ^ is given as 

e= (4( / i 1 - / i 2 ) 2 j(j+j) + i)^ 1/2 . (6) 

The upper and the lower signs in (4) correspond to the quasisinglet sg q and quasitriplet tr q 
states respectively. Note that g Si are the intrinsic spin magnetic moments of the constituent 
particles. According to the Dirac theory a free particle at rest has g Si = 2. In QED g Si 
is modified by the "anomaly", which to lowest order is given by the Schwinger correction 
g Si = 2 + a/n [7], where a is the fine-structure constant. 

In the case when m 2 » m\ our general result reduces to the previously known result 
[3-6], in which the orbital motion of the heavy particle is ignored, namely 

n _ „ F(F + l)+ 3l ( Jl + l)-/(/+l) 

9i 9h 2F(F + 1) ' {) 
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where 

9n ~ 1 + (9 S1 ~ 1) 2- (ji - 1} , (8) 

and 

„ _~ F(F+l)- 3l {ji + !) + /(/+!) . . 

2F(F + 1) • (9) 

Here J is the spin of the second particle, and F is the total angular momentum of the system. 
To facilitate the comparison of (4) with (7)- (9) we need to make the following replacement 
of quantum numbers: F — > J, J — > ji, l\ — £, S — > s±, I — > s 2 . 

In this paper we present results for the relativistic case beyond the formulae (2)-(4). 
Thus, the g-factors in (1) are written more generally as 

m = gf R + ^ E \ (io) 

where gf R are the nonrelativistic Lande factors defined by (2)-(4), and Ag REL is the rela- 
tivistic correction. In the next section we calculate /\gf EL to order O ((c/) 2 ) for all quan- 
tum states and arbitrary masses of particles. In most expressions we use natural units, i.e., 
h — c — 1, a — e 2 /An. The coupling constant is defined as a' = QiQ 2 /An. 

2. Variational wave equation and relativistic corrections to g- 
factors of two-fermion systems 

The relativistic wave equations for two-fermion systems in the absence of external fields 
were derived in [1] and [8] on the basis of a modified QED Lagrangian [9]- [10]. In this 
approach a simple Fock-space trial state of the form 

Atrial) = E/ rf 3 Pl^P2^ lS2 (pi,P2)^ lSl ^ 2S2 |0> , (11) 

S1S2 

is sufficient to obtain the HFS levels correct to fourth order in the coupling constant a'. Here 
& qm and -D q2S2 are creation operators for a free fermion of mass m\ and an (anti) fermion 
of mass m 2 respectively, and |0) is the trial vacuum state such that 6 qisi |0) = -D q2S2 |0) = 0. 
The F SlS , 2 are four-component adjustable functions. 



The variational principle S ( iptriai 



H -E 



Atrial / — 0, where H is the QED Hamiltonian 



is invoked to obtain a momentum-space wave equation for the amplitudes F SlS2 [1]: 

= S/ d 3 p 1 d 3 p 2 (u; pi +n p2 -E)F slS2 (p 1 ,p 2 )SF: iS2 (p u p 2 ) (12) 
d 3 pid 3 p 2 d 3 qid 3 q 2 



S1S2 



m 1 m 2 x - r 
~~ (2tt) 3 „ ti. J 



x ^ia 2 (qi,q2) (-i) M SlS2aia2 (pi,p 2 ,qi,q 2 ) SF* iS2 (pi,p 2 ), 
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where = p? + m\ and = + m^. The inter-particle interaction is represented by 
the generalized invariant .M-matrix, 

M SlS2ai(72 (pi, p 2 , qi, q 2 ) = M^ S2aia2 (pi, p 2 , qi, q 2 ) + M { ^ S2fTia2 (pi, p 2 , qi, q 2 ) + ••, (13) 

obtained as part of the derivation. It includes reducible and irreducible effects in all orders of 
the coupling constant a', and the sum contains all relevant Feynman diagrams. A discussion 
of the derivation and structure of the .M-matrix to one-loop level is provided in [1]. This 
equation allows one to obtain, in principle, all relativistic and QED corrections to the g- 
factor. 

The lowest-order QED corrections appear within the term and can be formally 

included in the intrinsic factor g Si , however we shall not do so in this work. In this paper we 
restrict our consideration to the first term Ai^ of the expansion (13), i.e., only tree level 
diagrams are included. The term Ai^ contains only relativistic corrections and it can be 
broken into two parts, namely 

M^ S2aia2 (pi, p 2 , qi, q 2 ) = M^ e S2aia2 (pi, p 2 , qi, q 2 ) + M e s f S2na2 (pi, p 2 , qi, q 2 ) , (14) 

where M° s ^ 2(7ia2 (pi, p 2 , qi, q 2 ) is the usual invariant matrix element, corresponding to the 
one-photon exchange Feynman diagram [1]. The element M. e ^ t S2ai(T2 represents the interaction 
with a given external classical field A^, 

(Pi,P2,qi,q 2 ) (15) 

= i (2nf 2 ( ^T^ A f(Pi ~ (Pi» *i) HQi) 7"« (qi, ax) 5 S 
\ +^^Af(q 2 - p 2 ) V (p 2 , a 2 ) (-iQ 2 ) YV (q 2 , s 2 ) 5 S 

Using the semi-relativistic expansion of the expression u(pi, Si) (qi, 0"i) up to order 1/c 3 
we obtain for the Ai — matrix 

>C 2 ^ 2 (Pi, P2, qi, q 2 ) = M^ 2 x X a2 (pi, p 2 , qi, q 2 ) + Mf^ ia2 (pi, p 2 , qi, q 2 ) , (16) 
where 

^i 1 1 £U(pi.p 2) qi ) q 2 ) (i7) 

(2tt) 3/2 ( £Mf( Pl - qi )pt (i (a, x ( Pl - qi )) + qi + P i), ^<W 3 (p 2 - q 2 ) 




2c \ + ^Af(q 2 - p 2 )xl 2 (i {<T 2 x (p 2 - q 2 )) + q 2 + p 2 ) . XsA^S 3 ( Pl - qi) 
is the non-relativistic contribution, and 

^£t 2 (pi.P2.qi.q 2 ) (is) 
' Q> a^Ht, n w ( (p?-q?)(qi- i ( (j i x qi)) \ m * n ) 

^A j ( Pl - qi )^ Sl ^ _ (p 2 _ q2) (pi + • ((Ji x pi)) ¥>ai<W (pa - q 2 ) 

, Q2 Aext (n _ W f ( (Pl ~ ql) (<J2 - i fa X q 2 )) \ 3 

V^A 5 (^-P2)x,2^_ (p 2_ q 2 )(p2 + i((72Xp2)) J X.W ( Pl - Ql ) 



2c 



V 
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(with ip\ = [1 0], y4 = [0 1], x\ = [0 1], xl = -f 1 °]> and J = 1,2) is the lowest-order 
relativistic correction to the non-relativistic term M.^ ext . 

For a stationary uniform magnetic field B = Bz the non-zero Fourier components of the 
vector potential are 

Afik) = £l& fe) » (*,) , -4f (k) = (y s %) . (19) 

The trial state (11) is taken to be an eigenstate of total linear momentum P, total angular 
momentum squared J 2 , its projection J 3 , parity V, and the Hamiltonian H, which corre- 
sponds to the hyperfine interaction [8]. In the rest frame, where the total linear momentum 
vanishes, the adjustable functions -F slS2 (pi, p 2 ) = F slS2 (p 1 )S (pi + p 2 ) can be specified for 
two categories of relations among the adjustable functions F slS2 (p): 

(i) The spin-mixed (quasi- singlet and quasi-triplet) states 

In this case we have £ = J, and the general solution under the condition of well-defined 
eigenvalues of P, J 2 , J 3 , and V can be expressed as [1], [8] 

F S1S , (P) = C£i msis2 /f ) { P )YT S ^ (P) + cS msiS2 /f \ P )Y^ ( P ) , (20) 

where m\\ = 1, m\2 = vri2\ = 0, m 22 = — 1. The Cj S ^ msiS2 = (imeSms \ Jmj) are the 
Clebsch-Gordan (CG) coefficients with total spin S, where S = (with index Si) for the 
singlet states and 5 = 1 (with index S* 2 ) for the triplet states respectively. Here fj Sl \p) 
and fj S2 \p) are radial functions to be determined. They represent the contributions of 
spin-singlet and spin-triplet states (the total spin S — 0, 1 is not conserved in general). 

(ii) The l-mixed triplet states 
These states occur for £i )2 = J =p 1. Their radial decomposition can be written as 

F S1S2 (P) = C IS2 ^WC 1!2 (P) + CjT S1 ' S2 f^K SlS2 (P) • (21) 

Again, the cf^™ 1 " 1 " 2 = (l-y^m^^Sms \ Jmj) are CG coefficients. For these states the 
system is characterized by J, mj, and V = (— 1) J . The orbital angular momentum £ = £ 12 
is not a good quantum number. Mixing of this type occurs only for states with principal 
quantum number n > 3. 

From the variational principle we obtain a system of coupled radial equations expressed 
in matrix form as 

(u, p + n p -E)W(p) = ^ I q2dq K(p, q )¥(q), (22) 

where a; 2 = p 2 + m\ and fi 2 = p 2 +m 2 , and q = |q|. Here F (p) is the two-component matrix 
of radial functions 



F(p) 









_ ff 2 \p) _ 


1 


_ Up) _ 



(23) 
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for spin-mixed and ^-mixed states respectively. The kernel of this equation is the 2x2 matrix 
[K] • • —fCij, which has the following form 



x {M° s i e S2ai(T2 ( P i, p 2 , qi, q 2 ) + A^SU (pi. Pa, qi> q 2 ) + MSToU (pi, P2, qi, q 2 )) • 

Here the are related to the CG coefficients by: C^f 2 = Cf^™' 1 ' 2 cfj?™' 1 ' 2 

and C^r 2 = cS msiS2 Cjm mCT1<T2 f o r the spin- and £-mixed states respectively. For the 



spin-mixed states we should take £i = £j = £. 

The solution of equation (22) with kernel (24) including only the first two terms with 
the M ope and M^ ext matrices was discussed in [2]. This solution describes the Zeeman 
splitting of the HFS energy levels in the non-relativistic limit. These energy levels are given 
by formulae (l)-(4). 

In order to obtain the Lande factors to order O (a' 2 ) we solve the radial equation (22) 
with the additional term M.^ ext in the kernel (24), which is evaluated perturbatively. The 
energy eigenvalues can be calculated from the matrix equation (22) as follows: 



In [2] we show that this system decouples for the spin-mixed states if the radial functions 



ff l] and /f l} are taken as = v / 0±0 IV J and /f l} = ^V(l T /2fj. Here 



fj = fi is a common radial function, the upper and lower signs correspond to sg q and 
tr q states respectively. The energy corrections for £-mixed states can also be calculated 
independently for £ = J — 1 and £ = J + 1 states with corresponding radial functions 
fe=j^i (see [11]). We evaluate (25) perturbatively using hydrogen-like radial functions (non- 
relativistic Schrodinger form fi = f^jm (p)) i n momentum space [4]). 

The calculations are straightforward, and yield the relativistic corrections to the g-factor 
for both particles of the system (the mass factors fa are defined by (5), indexes i and j are 
defined as in section 1). 

For the triplet states I = J — 1, 




(24) 




(25) 




(26) 



For the triplet states I = J + 1, 




(27) 
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For the spin-mixed states I = J, 

^ Ei =-<^.i^)(£) 2 > <»> 

where the upper and lower sign in (28) corresponds to sg q and tr q states respectively. For 
the spin-mixed states the parameter £ is given in Eq. (6). 

There are several works based on the Breit equation (e.g. [12]- [15]), where the relativistic 
corrections to the ^-factor of two-fermion bound states were considered. Some of the calcu- 
lations were obtained only for small mass ratio [12]; in others the orbital motion of the heavy 
particle was ignored [13]. These calculations did not take into account the mixed nature of 
the quantum states. We emphasize, that only mixed states diagonalize the Hamiltonian of 
HFS [2], [8]. Our result (26)- (28) is new and overcomes the above-mentioned shortcomings. 
For the relatively simple case of the ground state 15*1/2 (J — 1, V — —1) our results agree 
with the previous result of Hegstrom [14] and Grotch [15], namely Ag REL = — /^a /2 /3. Note 
that their definition of g e based upon AE e = ii B1 g e Bm s , Eq. (12) of [14], differs by a factor 
of 2 from our definition (1). 

We provide our result (26), (28) in numerical form for the lighter particle in atomic 
hydrogen (Table 1), muonium (Table 2), and muonic hydrogen (Table 3), for which a' = a. 
We consider only states with principal quantum number n = 1,2. For comparison, the 
nonrelativistic Lande factors g^ R ((2), (4)) with g Sl = 2 are also included in the tables. We 
used the following values for the mass ratios: m e /m p = 5.4461702 x 10~ 4 and m p /m^ = 
8.8802433. The fine-structure constant is taken as a = 7.2973523 x 10~ 3 . 

Table 1. g factor of the electron e~ in atomic hydrogen for n — 1,2 states. 
gi R - nonrelativistic Lande factor (2), (4), Agf £L -relativistic correction (26), (28) 



pe 


151/2(7=1) 


25i/ 2 (j=i) 


2-Pl/2(J=l) 


2-P3/2(j=i) 


2-P3/2(J=2) 


n NR 

9i 


1.0000000 


1.0000000 


0.3330513 


1.6661323 


0.9997278 


Ag KEL 


-0.0000177 


-0.0000044 


-0.0000133 


-0.0000133 


-0.0000053 



Table 2. g factor of the electron e in muonium for n — 1,2 states. 
gi R - nonrelativistic Lande factor (2), (4), Agf £L -relativistic correction (26), (28) 



H + e 


15l/2(J=l) 


25i/ 2 (j=i) 


2-Pl/2(J=l) 


2-P3/2(J=l) 


2-P3/2(J=2) 




1.0000000 


1.0000000 


0.3308504 


1.6619300 


0.9975935 


Ag REL 


-0.0000176 


-0.0000044 


-0.0000131 


-0.0000131 


-0.0000053 



Table 3. g factor of the antimuon /i in muonic hydrogen for n — 1,2 states, 
gj^-nonrelativistic Lande factor (2), (4), Agf -^-relativistic correction (26), (28) 



p + fx 


15l/ 2 (J=l) 


25i/ 2 (j=i) 


2-Pl/2(J=l) 


2-P3/2(J=l) 


2-P3/2(J=2) 




1.0000000 


1.0000000 


0.2880562 


1.5606937 


0.9496033 


Ag KUL 


-0.0000143 


-0.0000036 


-0.0000099 


-0.0000098 


-0.0000040 
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For the systems considered in Tables 1-3 the relativistic corrections for the heavier par- 
ticle Ag2 EL are negligible in comparison with Agf EL due to the small mass factor /ii. The 
relativistic corrections shown in the tables are small. They would be higher for the corre- 
sponding states in high-Z ions. For a realistic comparison of the (/-factor with experiment 
one would also need to calculate the QED corrections up to second order in a, that is the 
higher order matrix element Ai^ of Eq. (14) would have to be included [3]. 



3. Relativistic corrections to the g-factor in the one-body (Dirac) 
limit 

In this section we show the validity of our results for the g factor in the one-body limit. 
Note that the applicability of the formulae (2)-(4), (26)-(28) is restricted by the condition 
AE ext << AE HFS (i.e., a weak magnetic field B). To lowest order in a the HFS energy 
splitting for all states [8] (for m 2 >> m-y) is given as AE HFS fa a'V— . In the limit 
m 2 — > oo the HFS disappears, and in this case, the condition AE ext « AE HFS can not be 
satisfied for a nonzero magnetic field. 

In this case we need to go back to the original variational equation (12) and rewrite it in 
a form acceptable for the one-body limit. It is not difficult to show that with the trial state 

Atrial) =Y,f rf3 P^(P)t 1°) > ( 29 ) 
s J 

equation (12) reduces to the integral equation 
= W d 3 P K-E)F s (p)^(p)-^ j ^^F CT (q)(-i)^(p,q)5F;(p), 

(30) 

where E is the total one-body energy. In analogy to the two-body case, the matrix M S a (p, q) 
is made up of two parts (up to 0(a' 4 )). The first part corresponds to the one-photon 
exchange term M.°£ e S2(Tia2 (pi, P2, Qi, Q2) taken in the limit m 2 — > 00. This part describes 
the interaction of the particle with a static Coulomb potential A4°^ 2(Tl(T2 (pi, p 2 , qi, q 2 ) — > 
j^Couiomb jg-^ rj^g seconc j p ar ^ represents the interaction with an external magnetic field, 

namely M e s f (p, q) = M^J ext (p, q) + M.f} ext (p, q), which can be obtained from (17)-(18) 



M^ ext (p, q) = {2 l )3/2Q Af(p - q)v4 (i (a x (p - q)) + q + p) . Va 



2mc 

MfJ ext (p, q) = ( f£^ Af (p-q)yt (( p 2 _ q 2) (q _ , {a x q)) _ ( p 2 _ q 2) (p + , (c7 x p)) ) . ^ 

(32) 

Further calculations require a classification of the states. The trial state (29) is taken to be 
an eigenstate of total angular momentum squared j 2 and its projection j'3. These conditions 
can be satisfied if the adjustable two-component functions F s (p) are taken in the following 
form: 



(31) 
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For states £ = j-l/2 

Fi (p) = f)- h (p) C^Yp-? 2 (p) , F 2 (p) = /._, (p) C^'-'^ 2 (p) (33) 
For states i = j + l/2 

i^i (P) = (P) C 4 ^?i7 2 1/2 (P) » F - (P) = (P) C 4 "^?V2 1/2 (P) . (34) 

where Cj™* sr " s = (£mism s \ jrrij) are the CG coefficients for s = 1/2, m s = ±1/2. 

After substitution of these formulae into (30) and completion of the variational procedure 
we obtain the radial equation 

H - e) n (p) = -£^f -^L)c e (p, q ) n ( P ) , (35) 

where the kernel /Q (p, q) is expressed through the matrix M. and coefficients CJ^f. = 
q^CT*' namely 

K t (p, g) = -iJ2 c Jij (m%T* + A^«r (P, q) + A^SJT* (p, q)) 1^'* (P) n mCT (Q) , 

cr'cr 

(36) 

with a' = 1, 2, a = 1, 2, and mi j2 = rrij =F 1/2. 

In the absence of an external magnetic field, equation (35) represents the relativistic 
radial equation in momentum space for a bound one-body system. As was shown in [8], 
the solution of the two-body equation (22) reduces to the solution of the one-body equation 
(35). 

We now evaluate the contribution of the next two terms M.^ xt (p, q) and M^ xt (p, q) 
of equations (31) and (32). The relevant energy AE?^ e can be calculated perturbatively 
like in the two-body case (cf. Eq. (25)). We obtain 

AE% j>e = fiBgoBm,, (37) 

where 

9D = gL(e,j) + Ag$ SL {n,j). (38) 

Here ol (£) is the usual result for the anomalous Zeeman effect, being the Lande g factor, 
whose value is ([4]): 

tee.i) = f£r- (») 

The second term is the relativistic correction to the go factor 



Formulae (39) and (40) coincide with the result of the expansion (up to O (a' 2 )) of the general 
formula for the o-factor obtained by Margenau [16] on the basis of the Dirac equation with 
a' = Za. 
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4. Concluding remarks 

We considered the relativistic theory of the Zeeman splitting and the g-factor of the 
hyperfine structure in the two-fermion system on the basis of variational relativistic equations 
derived from quantum electrodynamics [1]. Relativistic corrections to the ^-factor beyond 
the non-relativistic formulae (2)-(4) (from Ref. [2]) were calculated up to order O (ct' 2 ) for all 
quantum states, and are given in Eqs. (26)-(28). The ^-factor corrections take into account 
the mixed nature of the states (spin-mixing and ^-mixing), and the orbital motion of the 
heavy particle. They are obtained for arbitrary mass ratio and are symmetrical with respect 
to the masses of the constituent particles. For the ground state our result reduces to the 
well-known formula obtained some time ago by Hegstrom [14] and Grotch [15]. We show 
that in the one-body case the solution of the variational relativistic equation reproduces the 
result for the ^-factor obtained from the Dirac equation. 

In our calculations we assumed that the trial state (11) is an eigenstate of the total 
linear momentum operator P. However this is only an approximation, because one can show 
that the operator P does not commute with the HFS Hamiltonian. To fix this problem we 
need to modify the trial state, or use an appropriate unitary transformation for the HFS 
Hamiltonian. The latter approach was discussed in the literature [12], [13], [17]. An analysis 
shows, that in our case the unitary transformation leads to the appearance of additional 
terms in the invariant .M-matrix. This is a technically difficult problem which we postpone 
for the future. 

We note that, in contrast to the Breit approach, the anomalous magnetic moment is not 
introduced from another calculation. As discussed in section 2 (below equation (13)) all 
QED effects are contained in the .M-matrix. The anomalous magnetic moment will appear 
naturally in our calculations if we include the next term A4.^ of the expansion of A4- matrix 
in (13). 

Concerning comparison with experiment we note that measurements so far appear to be 
restricted to states where the spin structure of the heavier particle can be ignored. This 
includes data for the 15*1/2 state in ions [18], [19]. The present work will be of practical 
importance when these measurements will be extended to all n = 2 (or higher-n) levels. 
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